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Abstract 

Corrections to the exact heavy-quark symmetry resuhs are expected to follow 
the l/rriQ mass effect of the heavy-quark. We show, by an explicit calculation, that 
there is something other than the mass effect that suppresses the breaking of the 
spin symmetry. 



1 Introduction 



The heavy-quark symmetry, which appears in the heavy-quark hmit, gives exact 
results for the decays of heavy hadrons |I| . Due to the heavy-quark symmetry all 
form factors in the heavy-to-heavy type of decays such as i? — > D^*^eh'e (-D*-*-* = D or 
D*) can be related, in the heavy-quark limit, to a single universal function called the 
Isgur-Wise function. The Isgur-Wise function is of nonperturbative origin and has 
been of great interest to both theoretical and experimental studies. In the heavy- 
quark symmetry limit, the decoupling of the heavy-quark spin with other light fields 
leads to symmetry relations among hadronic matrix elements. The corrections to 
the symmetry are expected to follow the l/mg mass effect of the heavy-quark. It 
has been unclear how well these would extrapolate to heavy-to-light quark decays, 
although presumably the charm quark might still be heavy enough. 

Many of these symmetries were anticipated in some versions of the constituent 
quark model but there has not been an estimate of how much of these were an artifact 
of the quark model, and of the choice of wave functions. In our recent paper using 
the relativistic quark model p[, we found that the breaking of the spin symmetry 
among hadronic form factors is small even for heavy-to-light quark decays. In this 
paper, we show explicitly, using the same model, that there is something other than 
the mass effect that suppresses the breaking of the spin symmetry. In fact, the quark 
model keeps the spin-symmetry rules remarkably well for a wide range of masses. 

We first recapitulate some aspects of spin symmetry for mesons in the heavy- 
quark limit. For a pseudoscalar meson P{Qq) with heavy constituent quark Q, the 
spin of Q decouples from all other light fields in P [Q. We can therefore construct 
the spin operator Sq for Q such that in the heavy tuq limit 

S§\Pm) = llVLiQq)) , (1) 

where VL{Qq) is the longitudinal component of a vector meson with the same quark 
content as P. In practice, the spin symmetry in Eq. (|1|) can be transformed into 
identities between the hadronic matrix elements, and thus some form factor relations, 
for H ^ P and H Vl, where H{hq) is a pseudoscalar meson with a heavy-quark 
h. Using the relation [| + A^] = (1/2) + V^) for the currents V;, = Q-f^h 

and A^ = (^7/^75^ of the transition h ^ Q, it can be shown that Eq. (|I]) leads to the 
following identity between the hadronic matrix elements for H ^ Vl and H ^ P, 

{Vl\A^ + A^\H) = {P\V° + V^\H) . (2) 

In \P) and \Vl) the spatial momentum of the quark Q is defined in the 2;-direction 
for the Q spinor to be an eigenstate of Sq. The spatial momenta of P and Vl 
should therefore also be defined in the z-direction such that the correction to the 
spin symmetry is of the order of A/mg, where A is the internal energy scale of P 
and Vl- 



2 Kinematics 



In this paper, we consider the breaking of the spin symmetry coming from a finite 
quark mass mg by directly calculating, in particular, the hadronic matrix elements 
in Eq. (Q). We use the relativistic quark model formulated in the infinite momentum 
frame (or equivalently, the light-front quark model) p|, ^, ^, |^ . We first define the 
ratio of the matrix elements in Eq. as 

, ^ {P{k')\V^ + V'\H{p)) 
' {VL{k)\A^ + A^\H{p)) ' 

so that 1 — p represents the percentage breaking of the spin symmetry. The ratio p is 
a function of momentum transfer such that (Vl|^° + A^\H) and {P\V^ + V^\H) are 
evaluated at the same g^. This allows us to evade the usual problem of kinematic 
discrepancies when we come to consider a number of different final states. Since the 
spatial momenta of P and Vl are in the ^-direction, we define the function p in a 
frame where the parametrization of the momenta in the 2;-direction is given by 



= {Eh; 0,0,f) , 
P = {Ev;0,0,k') , 

k"' = {Ep;0,0,k'') . (4) 

The vector and scalar masses my and mp are different for finite mq, so Vl and 
P will not carry the same momentum even though the initial state H has the same 
p^. We write the momenta k^ and k'^ in terms of the frame parameter p^ through 
the condition = {p — k)"^ = {p — k'Y- The general parametrization in includes 
the particular case of the infinite momentum frame in which k^ = k'^ = p^ = P 
where P ^ oo and = 0. It is important to note that the function p, when 
calculated in the infinite momentum frame, is defined at = only. This is also 
the point of maximum recoil, which is usually difficult to treat in a non-relativistic 
quark model, since a large amount of energy is given to the outgoing particle. In 
the infinite momentum or light-front frame, we have the following connection for p 
to the form-factors defined in Ref. pf; 

However, here we shall calculate p^^^^ directly from the matrix element. 

The mass-shell conditions for p, k, and k' give the following constraints on the 
momenta in M) 



where w = p-k/{mHmv) and w' = p-k'/{mHfnp). The ratio {Ep + k'^) / [Ey + k^) in 
(^) is therefore invariant for the frame defined in and is a function of through 
the relation = {p — k)"^ = [p — k'Y. 

It can be shown from the covariant expansion of the hadronic matrix elements 

that {P{k')\V'^ + V'^\H{p))/{Ep + k") and {VL{k)\A'^ + A'^lHip)) / {Ey + k') are 
invariant with respect to the frame defined in @. Using the kinematic constraint in 
(P), it is easy to see that the function p{q^) is an invariant quantity. The matching 

01 p^ va. H ^ Vl and H ^ P of p is the only choice that would lead to this 
invariance. In the definition of piq"^), there is an ambiguity coming from the fact 
that the ranges of q^ are usually quite different in H ^ Vl and H ^ P. We will 
therefore consider the value of p at = only. The hadronic matrix elements in 
(P) can be calculated reliably using the relativistic quark model formulated in the 
infinite momentum frame. So at = 0, we can write 

P(0) = , (7) 

where calculated in the infinite momentum frame. 



3 Symmetry Breaking 

A brief introduction to the relativistic quark model in the infinite momentum frame 
can be found in Refs. H, § . In the relativistic quark model, the wave function for 
the ground state meson M{Qq) is given by 

|M(k)) = / rfpQ E KrMPQ, ^)?lk - pQ, a)) , (8) 

where k = Pz is the spatial momentum of the meson M, pq = (p^-, xP) and 
Pg = k — pq = (— Pt, (1 — x)P) are those of the quarks Q and q, respectively, 
in the infinite momentum frame. Here, is the momentum wave function for 
the Qq bound state. It has the separable form into the spin and orbital parts as 
^aTcto- — ^iria ^1^1 1 where the expressions for R'{^'^^ and 0m can be found in Ref. 

In the relativistic quark model, it has been shown that p^*^^ is a function of the 
mass ratios r, s and /, so that p(0) = p^^^^{r, s, I) , where 

_ mg _ ^ _ A 

ruh ' rrih ' rrih ' 

The parameter A determines the internal energy scale of the meson and should be 
of the order of Aqcd- The dependence on / appears only in the momentum wave 
function and actually there could be a separate A for each of the mesons resulting in 
three further parameters. We take them all to be equal here. The kinematic region 
of interest for r, s, and / is such that < /, r, s < 1 . We find. 



(9) 



where 

Ji = / dx \ dyy(j)H(f)p — -rj- — ryj ^ (10) 

Jo Jo dn{l, sjaoir, s) 



and 



rl roo 

J2 = / dx dyy(f)H 
Jo Jo 



V- 



ao( l, g)ai(r, s)a2{r, s) + [ai(r, s) - a2(r, g) + ao(l, g) ] 
4(1, s)<ii(r, s)d2{r, s) 

(11) 

with the definitions 



Q;o(r, s) = xs + (1 — x)r , ai(r, s) = r + xMo(r, s) , a2(?", s) = s + (1 — x)Mo(r, s) , 

^o(^s) = W + - , 

V X 1 — X 

cio(r, s) = \/al{r, s) + y^ , di{r, s) = \J aj{r, s) + y^ , diir, s) = \/ al{r, s) + y^ . 

We may write the orbital wave functions (pn and 4>py in terms of a Gaussian 
function 0(r, s, I) such that = 0(1, s, /) and (ppv = s, I). The expression for 
0(r,5,Oisgivenby 

<P{r,s,l) = N^exp^-^iy' + z')/l'^ , (12) 

where is a normahzation factor that is canceled out in p^^'^^ , and 

In,./ n (r^-s^) 
z = (x- -)Mo(r, s) - ; — { . 

In Ref. 0, it has been pointed out that the scaling behavior of the meson decay 
constant Jm in the heavy-quark limit imposes a constraint on the orbital wave 
function. The Gaussian wave function in ([I^ ) is shown to satisfy the scaling law 
1/ y/rriQ of Jm in the heavy rriQ limit. 

In the numerical analysis of p(0), it is convenient to set s = Z and vary s and 
r within the kinematic region of < s,r < 1. The spectator quark q is thus 
considered to be a light quark with rriq = A ~ Aqcd- In case of a heavy-quark 
decaying to a heavy or to a light quark, the corresponding regions for s and r are 
such that s is small and r varies between 1 and 0. By comparison, for a light quark 
decay to another light quark, we look at the region where s and r could both be close 
to 1. When s = r = 1, the symmetry breaking is calculated to be 1 — p(0) = —0.33, 
using the Gaussian wave function in (0). Thus the heavy-quark symmetry does 
not apply in that C3iSG, clS expected. 

In Fig. (|l|), we show the plot of 1 — p(0) for a very heavy-quark decaying to a 
heavy or to a light quark so that the the mass ratios s and / are very small (here they 
are taken to be s = Z = 0.001). The variation of r is within the range < r < 1. 
From the Gaussian form of the wave function in Eq. ([T2|) one expects {y"^) = and 



from the expression for p^^^^ in Eq. (y), we expect 1 — p(0) — * since (y^) for 
a very heavy decaying quark. 

In Fig. (|I|), we show that the breaking 1 — p(0) is less than 1% in the heavy- 
quark hmit (or small s = 0.001). From the figure, the mass effect of rriQ can be seen 
clearly as the breaking of the spin symmetry gradually increases from large r, or 
heavy-to-heavy decays, towards the smaller r, or heavy-to-light region. However, 
even at very small r, corresponding to a decay b ^ s {r = 0.1) or b ^ u {r = 0.06), 
it is remarkable that the symmetry breaking is less than 0.6%. 

In Fig. (|^), we show the plot of 1 — p(0) with the physical spectator quark 
mass ratios for heavy-quark decays as s = rriu/mb = 0.06, s = ms/mi, = 0.1, 
s = rriu/iTLc = 0.2, and s = ms/mc = 0.3. In the figure, the breaking of the spin 
symmetry is shown to be less than about 10% for heavy b and c quark decays. In 
the strict heavy mass limit of s ^ as indicated in Fig. (|I]) , the mass effect of 
rriQ is clearly seen as the suppression of the symmetry breaking with increasing r. 
Notice, however, that for finite spectator mass ratios, the function 1 — p(0) passes 
through zero at a recoil mass r = mg/m/j below its heaviest limit r = 1. 

In the case of b decays and the decay of charm to non-strange quarks, the mass 
effect of mgis still seen in the region where 1 — p(0) is positive for a large range of 
r. In the region where 1 — p(0) is negative, the ttiq suppression no longer follows 
as the size of 1 — p(0) increases with r = mQ/mh. This suggests that there are 
kinematic factors other than the mass effect of uiq that govern the size of the 
symmetry breaking 1 — p(0), when the decaying quark has finite mass. As shown 
in the figure, the zero of 1 — p(0) is at smaller r when the ratio s is larger. The 
zero for s — appears at r — > 1 and decreases with increasing s. For s > 0.251, 
1 — p(0) is negative for all r. The mass effect of mq is therefore less pronounced as s 
gets larger and the kinematic effects dominate. A different type of behavior enters 
for the decay involving a charm quark and a strange spectator where the deviation 
from the symmetry limit is rather constant and about 4%. 

In fig. d^), we show the corresponding values of s and r for which 1 — p(0) = 0. 
The zero is shown to lie within the region where s is small {s < 0.251) and the 
decaying quark is heavy. As shown in the figure, the zero of 1 — p(0) appears at 
smaller r as the mass of the decaying quark becomes less heavy (larger s). Also 
shown in the figure are the plots for which the spin-symmetry breaking is about 
10% that is |1 — p(0)| = 0.01. It can been seen that a large portion of the possible 
phase space of r and s is within the region where |1 — p(0)| < 0.01. 

We have also obtained a similar result using the harmonic oscillator wave function 

instead of the Gaussian function. The result is therefore not an artifact of a 
particular momentum wave function. The quantity p(0) has the following physical 
meaning: 

I (n\\2 _ i^H - ^vf dV{H Plu)/dq\2=Q 
'^^ ^' {w?H-mlfdV{H^VLlv)/dq^\^ " ^ ' 



,1 



This allows a test of these results to be made by considering the g spectrum for the 
semileptonic decays H{hq) — > P,VL{Qq)- The size of p(0) for particular values of 



r and s can now be measured. Repeating this for the different semileptonic decay 
channels of H, the dependence of p(0) with r and s can also be determined. 
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Figure 1: The plot of 1 — p(0) using the Gaussian orbital wave function. The plot is 
for heavy-to-heavy and heavy-to-light decays with a quark mass ratio s = I = 0.001 
and where r varies within the range of < r < 1 . (r refers to the ratio r = ttiq / rrih 
in the quark decay h ^ Q). 
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